meta %>%
filter(su_blkgp == 1) %>%
select(varname, about) %>% as.list()
## $varname
## [1] "Blkgr" "consproviders" "busproviders" "avgMaxAdDown"
## [5] "avgMaxAdUp" "p_under25.3mbps"
##
## $about
## [1] "12-digit block group code"
## [2] "The number of consumer internet providers in the su"
## [3] "The number of business internet providers in the su"
## [4] "The average maximum advertised download speed by each broadband provider in the su"
## [5] "The average maxmum advertised upload speed by each broadband provider in the su"
## [6] "The proportion of census blocks within the block group/tract that have internet speeds under 25/3 Mbps, and therefore don't meet the FCC benchmark for \"advanced telecommunications capability\""
glimpse(cvldat)
## Rows: 155
## Columns: 6
## $ Blkgr <dbl> 510030101001, 510030101002, 510030101003, 510030102011…
## $ consproviders <int> 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 7, 7, 10, 8, 8, 8, 7, 7,…
## $ busproviders <int> 8, 8, 10, 9, 10, 8, 8, 10, 8, 13, 10, 8, 12, 9, 12, 10…
## $ avgMaxAdDown <dbl> 32.07490, 28.82663, 46.93761, 66.48129, 72.78802, 101.…
## $ avgMaxAdUp <dbl> 2.024672, 1.738201, 4.158629, 4.804956, 3.529576, 4.13…
## $ p_under25.3mbps <dbl> 0.25882353, 0.21052632, 0.18095238, 0.29213483, 0.1470…
cvldat %>% select(consproviders, busproviders, avgMaxAdDown, avgMaxAdUp, p_under25.3mbps) %>%
select(where(~is.numeric(.x))) %>%
as.data.frame() %>%
stargazer(., type = "text", title = "Summary Statistics", digits = 2,
summary.stat = c("mean", "sd", "min", "median", "max"))
##
## Summary Statistics
## ===================================================
## Statistic Mean St. Dev. Min Median Max
## ---------------------------------------------------
## consproviders 7.43 1.00 5 7 10
## busproviders 9.10 1.76 6 9 15
## avgMaxAdDown 112.59 59.70 17.99 105.60 239.96
## avgMaxAdUp 27.04 33.63 1.56 8.41 173.55
## p_under25.3mbps 0.11 0.16 0.00 0.05 1.00
## ---------------------------------------------------
cvldat %>% select(c(Blkgr, consproviders, busproviders, avgMaxAdDown, avgMaxAdUp, p_under25.3mbps)) %>%
pivot_longer(-Blkgr, names_to = "measure", values_to = "value") %>%
ggplot(aes(x = value, fill = measure)) +
scale_fill_viridis(option = "plasma", discrete = TRUE, guide = FALSE) +
geom_histogram() +
facet_wrap(~measure, scales = "free")
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
meta %>%
filter(varname %in% c("consproviders", "busproviders", "avgMaxAdDown", "avgMaxAdUp", "p_under25.3mbps")) %>%
mutate(label = paste0(varname, ": ", about)) %>%
select(label) %>%
as.list()
$label [1] "consproviders: The number of consumer internet providers in the su"
[2] "busproviders: The number of business internet providers in the su"
[3] "avgMaxAdDown: The average maximum advertised download speed by each broadband provider in the su"
[4] "avgMaxAdUp: The average maxmum advertised upload speed by each broadband provider in the su"
[5] "p_under25.3mbps: The proportion of census blocks within the block group/tract that have internet speeds under 25/3 Mbps, and therefore don't meet the FCC benchmark for "advanced telecommunications capability""
pal <- colorNumeric("plasma", reverse = TRUE, domain = mapdat$consproviders)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(consproviders),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T
),
popup = paste0("GEOID: ", mapdat$Blkgr, "<br>",
"Number of internet providers: ", mapdat$consproviders)
) %>%
addLegend("bottomright", pal = pal, values = mapdat$consproviders,
title = "Number of internet <br> providers", opacity = 0.7)
pal <- colorNumeric("plasma", reverse = TRUE, domain = mapdat$avgMaxAdDown)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(avgMaxAdDown),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T
),
popup = paste0("GEOID: ", mapdat$Blkgr, "<br>",
"Average max advertised <br> download speeds: ", round(mapdat$avgMaxAdDown, 2))
) %>%
addLegend("bottomright", pal = pal, values = mapdat$avgMaxAdDown,
title = "Average maximum <br>advertised <br>download speeds", opacity = 0.7)
pal <- colorNumeric("plasma", reverse = TRUE, domain = mapdat$avgMaxAdUp)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(avgMaxAdUp),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T
),
popup = paste0("GEOID: ", mapdat$Blkgr, "<br>",
"Average max <br> advertised upload speeds: ", round(mapdat$avgMaxAdUp, 2))
) %>%
addLegend("bottomright", pal = pal, values = mapdat$avgMaxAdUp,
title = "Average maximum <br>advertised <br> upload speeds", opacity = 0.7)
pal <- colorNumeric("plasma", reverse = TRUE, domain = (mapdat$p_under25.3mbps *100))
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal((mapdat$p_under25.3mbps *100)),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T
),
popup = paste0("GEOID: ", mapdat$Blkgr, "<br>",
"Percent of blocks <br>without advanced <br>telecommunications capability: ", round((mapdat$p_under25.3mbps *100), 2))
) %>%
addLegend("bottomright", pal = pal, values = (mapdat$p_under25.3mbps *100),
title = "Percent of blocks <br>without advanced <br>telecommunications capability", opacity = 0.7)